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Abstract. We present high resolution optical spectroscopy of three candidate members of the Taurus- Auriga star 
forming region. Based on the spectral type, the strength, profile and width of the Hq line, the detected lithium 
at 6708 A, the location of these objects in a H-R diagram and the comparison with similar objects belonging to 
young stellar associations, we determine that they are bona fide members of the SFR, with about ~3 Myr, have 
masses at or below the substellar limit and, at least in one case, there is active accretion from a circum(sub)stellar 
disk. This result suggests that high mass brown dwarfs go through a Classical TTauri phase and form like stars, 
from colapse and fragmentation of a molecular cloud. 
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1. Introduction 

Brown dwarfs, objects unable to fuse hydrogen in an sta- 
ble manner (i.e., with masses below at about 0.072 M©, 
Baraffe et al. 1998), pose an important problem to the 
theory of stellar formation. Several formation mechanisms 
have been proposed, including formation like a star (from 
collapse and fragmentation of a molecular cloud, Padoan 
P., & Nordlund 2004) to a planet-like process (from a 
circumstellar disk) or as stellar "embryos", ejected from 
multiple systems before they are able to accrete enough 
matter (Reipurth & Clarke 2001; Bate et al. 2002). These 
proposed mechanisms have different implication in their 
formation, their evolution and their properties. In partic- 
ular, if brown dwarfs are created like stars, it seems that 
they should go through an phase of active accretion from 
a circum(sub)stellar disk, such as low mass stars do (Shu 
et al. 1987). Actually, since the last couple of years, dif- 
ferent groups have presented indirect and direct evidences 
of active accretion in a handful of young brown dwarfs 
belonging to several nearby star forming regions (SFR), 
open clusters and moving groups (Fernandez & Comeron 
2001; Muench et al. 2001; Natta & Testi 2001; Natta et al. 
2001, 2002; Testi et al. 2002; Jayawardhana et al. 2002ab, 
2003ab; MuzeroUe et al. 2003; Barrado y Navascues et 
al. 2002, 2003, 2004; Barrado y Navascues & Martin 
2003; Mohanty & Basri 2003; Mohanty, Jayawardhana & 
Barrado y Navascues 2003; Comeron et al. 2003) 

In this paper, we enlarge the sample by collecting high 
resolution spectroscopy of three proposed members of the 
Taurus- Auriga SFR, located at 140 pc and about 1-3 Myr 
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(for an update on Taurus, see Luhman et al. 2003 and 
references therein). 



2. Observations 

Our Taurus targets were selected from Bricefio et al. 
(2002). In that paper, a new sample of nine mid- to late- 
M candidate members were presented, based on Ic and 
2;' photometry (plus 2MASS near IR data) and low res- 
olution spectroscopy. We have collected high resolution 
spectra of a third of them -the high mass brown dwarfs, 
whose spectral type is between M5.75 and M7.5- with the 
Magellan I 6.5m telescope and the MIKE echelle spectro- 
graph on 2002, Dec ll-14th. Although three objects are 
a small number, they represent a significant fraction of 
the high mass brown dwarfs belonging to this association 
discovered so far (about 12 in the quoted spectral range). 
Therefore, they can provide some hints about their forma- 
tion mechanism. 

Additional details about the observations can be found 
in Barrado y Navascues, Stauffer & Jayawardhana (2004), 
where we analyze a sample of very low mass stars and 
brown dwarfs belonging to the ~5 Myr cluster associated 
to the A Orionis star. In order to optimize the signal-to- 
noise, we binned the data during the read-out to two by 
two pixels in the spatial and the spectral direction, respec- 
tively, yielding a spectral resolution of R=25,000 (^0.25 
A), with a 0.75 arcsec slit. The spectral range of our spec- 
tra is 4500-7250 A. 
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Fig. 1. Spectra around HeI6678 A and LiI6707.8 A. 

3. Analysis and Discussion 

3.1. Spectral types 

Spectral types were derived by comparing with several 
spectral templates, by using the order around 7050 A, 
corresponding to a TiO band, following Mohanty et al. 
(2004), since this range is very sensitive to effective tem- 
perature. Errors can be estimated as half a subclass. Our 
values are very close to those obtained by Briccfio et al. 
(2002) from low resolution spectra (see Table 1). 

3.2. On the lithium abundance, age, and mass 

We have detected lithium 6707.8 A in two out of the three 
targets (KPNO-Tau-05 and KPNO-Tau-08, Figure 1). The 
spectrum of KPNO-Tau-03 has worse quality. With some 
caveats, the visual inspection indicates that this feature is 
also present. Note that Briceho et al. (2002) states that 
they detected lithium with their low resolution spectrum 
(at higher signal- to- noise ratio). This element is easily 
destroyed in the stellar interior, being its surface abun- 
dance dependent on mass and age (as well as other sec- 
ond order parameters). In fact, brown dwarfs which are 
more massive than about 0.060 Mq do deplete lithium 
in a time scale of few tens of million years (D'Antona & 
Mazzitelh 1994; Burrows et al. 1997; Baraffe et al. 1998). 
The Taurus-Auriga complex has an age between 1 and 
3 Myr. Since our three candidate members have spectral 
types between M6 and M7, they should have masses equal 
or larger than 0.06 Mq, according to Baraffe et al. (1998). 
Therefore, the detection of this alkali clearly indicates 
that these objects are in the pre-Main Sequence (PMS). 
Moreover, from the statistical point of view, the likelihood 
of having three late-M, PMS interloper whose spectral and 
photometric properties coincide with the Taurus sequence 
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Fig. 2. Lithium equivalent width versus the spectral type. 
The solid line corresponds to the upper envelope of the 
values measured in young open clusters. The long-dashed 
line delimits the areas for weak-line and post-T Tauri stars 
(adapted from Martin 1997 and Martin & Magazzii 1999). 
Short-dashed and dotted lines correspond to the cosmic 
abundances -A(Li)=3.1~ from gravities of Logg=4.5 and 
4.0, respectively (curves of growth from Zapatero Osorio 
et al. 2002). All Pleiades and IC2391 members with mea- 
sured lithium equivalent width are shown as open circles 
and triangles -upper limits-. Note the lithium depletion 
boundary at ~M5.5. Sigma Orionis low mass stars and 
brown dwarfs appear as solid squares (Zapatero Osorio et 
al. 2002). The big asterisks represent KPNO-Tau-3, 5 and 
8. 

is negligible. Therefore, we have to conclude that they, 
indeed, belong to the association. 

Figure 2 displays values of the lithium equivalent width 
(W) measured for our targets (large asterisks), pre- main 
sequence mcmbcirs of the ^5 Myr Sigma Orionis clus- 
ter (Zapatero Osorio et al. 2002, solid squares), and stel- 
lar and substellar members of IC2391 -53 Myr- and the 
Pleiades 125 Myr- from Barrado y Navascues et al. 
(1999, 2004), Soderblom et al. (1993), Garcfa-Lopez et 
al. (1994), Jones et al. 1996, Stauffer et al. (1998) and 
Jeffries et al. (1999). The solid line describe the maxi- 
mum W(Li) measured in cluster stars. The long-dashed 
curve is an update version of the criterion used by Martin 
(1997) and Martin & Magazzii (1999) to distinguish be- 
tween Weak-line and post- TTauri stars. Finally, the dot- 
ted and short dashed lines correspond to the Zapatero 
Osorio et al. (2002) curves of growth for Logg=4.0 and 
4.5 and an abundance of A(Li)=:3.1 (i.e., cosmic abun- 
dance). We have labeled in the plot different sections. The 
W(Li) for KPNO-Tau-5 is lower than the measured values 
in the other two Taurus objects. This fact might be due 
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to the presence of optical veiling, about r(6700)~0.2, al- 
though we do not think this is the case. On one hand, no 
other sign of accretion has been found in this object (see 
next section). Moreover, this equivalent width is compat- 
ible with the dispersion found in Sigma Orionis and the 
-undcplctcd- values characteristic of cluster brown dwarfs. 
For the other two objects, no veiling seems to be present 
either. In the lithium equivalent width dispersion is real 
both in Taurus and Sigma Orionis cluster (the same effect 
is present in the Lambda Orionis association, Barrado y 
Navascues, Stauffer & Jayawardhana 2004), it might im- 
ply diferences in the surface abundances due to additional 
mixing mechanisms in the stellar and substellar interior, 
different to pure convection (such as those proposed in so- 
lar mass stars), differences in the structure due to differ- 
ent stellar and substellar parameters (such as fast rotation 
in the case of pre-Main sequence stars, Martm & Claret 
1996), or to induced effects related to activity and/or ro- 
tation over the spectral feature itself or the surrounding 
continuum (see the discusion in Barrado y Navascues et 
al. 2001 and references therein). 

As a summary, based on these data, mainly from the 
detection of lithium and accretion when present (see next 
section), we can conclude that these three objects are 
young and, indeed, belong to the Taurus region. In the 
case of KPNO-Tau-5, due to its spectral type (M7) its 
substellar nature is well established. The other two are lo- 
cated at the substellar borderline and their nature is not 
so firmly established due to uncertainties in the models 
and the spectral type determination. 

Once membership to the stellar association has been 
proved, we have derived the bolometric magnitudes from 
Ic and Ks, the distance modulus (m — Af)o=5.731 (140 
pc, Kenyon, Dobrzycka & Hartmann 1994), the redden- 
ing derived by Bricefio et al. (2002) and the bolomet- 
ric corrections by Comeron et al. (2000) and Tinney et 
al. (1993) for these two bands {Ic and Ks, respectively). 
Masses were computed based on models by Baraffe et al. 

(1998) . Effective temperatures were obtained using sev- 
eral scales, namely Luliman (1999) for intermediate grav- 
ity and Leggctt (2000, 2001). All the measured and derived 
values are listed in Table 1. 

Figure 3 displays a HR diagram. The isochrones and 
evolutionary tracks -solid and dashed lines, respectively- 
are from Baraffe et al. (1998). Asterisks and open stars 
correspond to bolometric luminosities obtained from Ic 
and A'.s. respectively. For this particular diagram, we 
made use of the effective temperature scale by Luhman 

(1999) . Other temperature scales, such as that from 
Leggett et al. (2001), would shift the location of these 
three objects (two of them are almost on top of each 
other, the small shift in Teff is arbitrary) to the right hand- 
side, making them younger and less massive. Note, how- 
ever, that Luhman' scale has been tunned specifically for 
the model we are using here. In any case, regardless the 
election of models (alternative models are, for instance, 
Burrows et al. 1997; D'Antona & Mazzitelli 1994, 1997, 
1998; Chabrier et al. 2000; Baraffe et al. 2002), bolometric 
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Fig. 3. HR diagram of KPNO-Tau-3 and 8 (on top of each 
other) and KPNO-Tau-5. The values represented by aster- 
isks and open stars were derived from Ic and Ks magni- 
tudes, respectively. 

corrections and effective temperature scale, these members 
are at or below the substellar frontier and have an age of 
about 3 Myr. 

3.3. Near infrared photometry, Ha emission and 
accretion 

None of our three targets seems to have neither the pres- 
ence of forbidden lines, characteristics of outflows, nor 
near infrared infrared excesses, coming from an accretion 
disk. Figure 4 compares the colors {Ic — J) and {H — Ks), 
and includes the loci for dwarfs -BesscU & Brett 1988; 
Kirkpatrick et al. 2000; Leggett et al. 2001)- and Classical 
TTauri stars (Meyer et al. 1997, Barrado y Navascues et 
al. 2003) and Classical TTauri stars from Orion (Herbig & 
Bell 1988). As can be seen, the Taurus objects have pho- 
tometric properties -data from 2MASS- similar to those 
MS stars of similar spectral type and, as stated in the 
previous paragraph, no near IR excess is seen. This fact, 
by itself, cannot prove or disprove the presence of a cir- 
cum(sub)stellar disk, since a hole can be present or the 
disk temperature can be very cold. Moreover, the IR ex- 
cess might depend on the orientation of the disk. See, 
for example, the case of LS-RCrA 1, M6.5 brown dwarf 
belonging to RCrA dark cloud (Barrado y Navascues, 
Mohanty & Jayawardhana 2004). Additional observations 
and longer wavelengths -more sensitive to cooler disks- 
both from ground-based telescopes and space-borne in- 
struments such as those in Spitzer Space Telescope, can 
help to shed some light in this issue. 

We have detected and measured Ha in emission for 
these three objects. The profiles of this feature are dis- 
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Fig. 4. Color-color diagrams. The position of the Taurus 
BD candidates are indicated with large asterisks. Crosses 
indicate the position of classical TTauri stars belonging to 
Orion stellar population (Herbig & Bell 1988). The thick- 
solid and dashed lines correspond to the locii of the main 
sequence stars (from Bessell & Brett 1988; Kirkpatrick et 
al. 2000; Leggctt et al. 2001) and CTT stars (Meyer et al. 
1997, Barrado y Navascues et al. 2003), respectively. 
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Fig. 5. Ha profiles. The instrumental profile is included 
as a dashed line (top spectrum). 



played in Figure 5. Note the possible asymmetry in 
KPNO-Tau-8, the double peak in KPNO-Tau-5 and 
KPNO-Tau-3, and the width of the line (310 km/s), 
typical of accreting brown dwarfs (White & Basri 2003; 
Jayawardhana et al. 2003). In fact, this last object is 



Fig. 6. Ha equivalent widths for members of the Taiirus 
SFR. Solid and open circles correspond to Classical and 
TTauri stars. Objects with no classification are shown as 
crosses. Large open circles or square represent objects 
with mid-IR excesses and forbidden lines, respectively. 
The three objects studied here are displayed as large as- 
terisks. The dotted, bold curve is the saturation criterion, 
whereas two previously proposed criteria (5 and 20 A) to 
separate CTTS and WTTS are included as long-dashed, 
thin horizontal segments. The vertical dotted segment de- 
notes the location of the substellar frontier. 



above the criterion defined by Barrado y Navascues & 
Martin (2003) which discriminate between accreting and 
non accreting objects, as Figure 6 clearly indicates. We 
note, however, that this criterion depends on low reso- 
lution spectroscopy, which normally yields larger equiv- 
alent widths compared with higher resolution data. In 
this diagram, we have included Classical and Weak-line 
TTauri stars belonging to Taurus as solid and open circles. 
Members without classification are included as crosses. 
Our three targets appear as large asterisks. Overlapping 
squares and big circles denote those members with forbid- 
den emission lines and near-infrared excesses, respectively. 
This large W(Ha) agree with the fact that we also detect 
HeI6678 A, another accretor indicator (see Figure 1). 

Active accretion in at least one Taurus member whose 
mass is close or below the substellar limit is important 
for several reasons. First, although the sample studied 
in this paper is very small, it suggests that a significant 
fraction of the very low mass stars and high mass brown 
dwarfs might harbor an accretion disk. Second, the new 
data amass additional evidence for accretion in the sub- 
stellar domain. To the best of our knowledge, there are 
only three papers dealing with high resolution spectra in 
Taurus brown dwarfs (White & Basri 2003; Muzerolle et 
al. 2003; and Jayawardhana et al. 2003). The first work 
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found three accretors in a sample of ten very low mass 
stars and brown dwarfs, although the less massive, a M6.5 
(GM Tau), has a mass just above the substcUar limit. 
The second paper includes four objects whose spectral 
type is M6 or M7. None of them seem to undergo ac- 
tive accretion based on the width of Ha, although the 
authors claim that one of them, namely MHO-5, is ac- 
creting based on the detection of forbidden lines of oxy- 
gen at 6300 and 6363 A and Call IRT. Moreover, three 
out of the four were observed at moderate resolution, in- 
cluding MHO-5 (R~8000). Regarding the later study, the 
four M7-M7.75 brown dwarfs discussed there, with masses 
down to 0.05 Mq (again, using Baraffe et al. 1998 mod- 
els), do not show accretion either. MuzeroUe et al. (2003) 
also analyzed seven Taurus members whose spectral types 
are M4.75-M5.75 (just above or at the substellar border- 
line). Three out of these seven are accreting, based on the 
width of Ha and other spectral indicators. Therefore, up 
to date, our study presents the only brown dwarf belong- 
ing to Taurus (kpno-tau-03), at the substellar limit, which 
has been proved to be accreting. 

In conjunction with the studies quoted in the previ- 
ous paragraph, our results indicate that about 10% of the 
Taurus brown dwarfs (one out of 11) is actively accreting. 
This fraction is much smaller than the accretion occur- 
rence in low mass stars in the association or the estimate 
for the substellar domain (about 50%) based purely on 
the strength of Ha measured in low resolution spectra 
(Barrado y Navascues & Martin 2003). The statistical cri- 
terion defined in this last work is based on the saturation 
of the activity, as measured in several young open clusters 
(namely IC2391, Alpha Per and the Pleiades, with ages 
ranging from about 50 to 125 Myr). The discrepancy in 
the fraction of accreting brown dwarfs might imply that 
there is an additional source of flux in Ha line or that 
the accreting criterion based on the Ha width at 10% of 
the maximum intensity, as defined by Jayawardhana et al. 
(2003), i.e. 200 km/s, is too restrictive. 

The detection of accretion in substellar objects in- 
dicates that they undergo a phase similar to Classical 
TTauri stars (for a review, see Bertout 1989 or 
Appenzeller & Mundt 1989). A handful of other brown 
dwarfs belonging to other young stellar association, with 
ages ranging from 1 to 10 Myr, have been observed with 
high resolution spectroscopy and the presence of accre- 
tion confirmed. In at least one, LS-RCrA 1 (Fernandez 
& Comeron 2001; Barrado y Navascues, Mohanty & 
Jayawardhana 2004), outflows have been seen, by means 
of the detection of intense, narrow forbidden lines. In a 
previous paper we have named them Classical TTauri sub- 
stellar analogs (CTTSA). Therefore, we can conclude, at 
least in the case of high mass brown dwarfs (Mas~0.072- 
0.04 Mq), that they are fornic;(l as low mass stars, by frag- 
mentation and collapse of the original molecular cloud. Of 
course, additional studies are needed to confirm this pre- 
liminary conclusion, in particular high resolution imag- 
ing in the near and mid- infrared, as well as in the opti- 
cal. Narrow band imaging might show whether these ob- 



jects present outflows similar to those observed in some 
Classical TTauri stars. 

4. Summary 

We have collected high resolution spectroscopy for the 
members of the Taurus-Auriga SFR. They have M6-M7 
spectral types with moderate to intense Ha emission, with 
profiles that, at least in one case, can be classified as typ- 
ical of an accretor. We detect lithium too, which indicate 
that they are in the pre-main sequence (more strictly, 
the equivalent width imposes a maximum age of about 
30 Myr, based on Baraffe et al. 1998). Moreover, at least 
one, possibly two of them, are accreting material, as shown 
by the Ha profile, providing an additional, stronger con- 
strain to the age, about 10 Myr. Therefore, a minimum of 
alO% (one out of eleven) of the high mass brown dwarfs 
known in the cluster show accretion (i.e., circumsubstelalr 
disks). Note, however, that these objects do not have near- 
infrared excess and do not show forbidden lines, although 
HeI6678 A, another signpost of accretion, appears in the 
spectrum of one of them, possibly in the other two. Since 
the likelihood of having a stellar association younger than 
30 Myr in the same line-of-sight as Taurus, we can safely 
conclude that these objects belong to this SFR. Based on 
the new spectral types and magnitudes in the optical and 
infrared, we derived masses around the substellar limits 
and an age close to 3 Myr. Finally, the detection of accre- 
tion in the substellar domain provides another clue about 
the formation mechanism of high mass brown dwarfs, sug- 
gesting that they form as stars, by fragmentation and col- 
lapse of a molecular cloud. 
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